Single-photon excitation of a coherent state: catching the elementary step of 

stimulated light emission 
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When a single quantum of electromagnetic field excitation is added to the same spatio-temporal 
mode of a coherent state, a new field state is generated that exhibits intermediate properties between 
those of the two parents. Such a single-photon-added coherent state is obtained by the action 
of the photon creation operator on a coherent state and can thus be regarded as the result of 
the most elementary excitation process of a classical light field. Here we present and describe 
in depth the experimental realization of such states and their complete analysis by means of a 
novel ultrafast, time-domain, quantum homodyne tomography technique clearly revealing their non- 
classical character. 

PACS numbers: 42.50.Dv, 03.65.Wj 



INTRODUCTION 

A coherent state \a), the eigenstate of the photon anni- 
hilation operator a\a) = a ja), is the closest analogue to 
a classical light field and exhibits a Poisson photon num- 
ber distribution with an average photon number of jap 
and a width of \a\. Coherent states posses well defined 
amplitude and phase, whose uncertainties are the mini- 
mum permitted by the Heisenberg uncertainty principle. 
On the contrary, a Fock state, the eigenstate of the pho- 



ton number operator (V a\n 



contains a perfectly 



defined number of quanta of field excitation and is strictly 
quantum-mechanical, with no classical analog. Moreover, 
being the intensity of a Fock state defined without un- 
certainty, its phase is completely undefined. 

Single-photon Fock states have been recently gener- 
ated by means of conditional preparation techniques 
and homodyne tomography has been used to completely 
characterize them with the reconstruction of an asso- 
ciated Wigner function clearly exhibiting negative val- 
ues U 0- Displaced Fock states, obtained by mixing 
a coherent state with a single photon upon a highly re- 
flecting beam splitter as in Ref . Q , have also been inves- 
tigated with tomographic techniques which have shown 
the non-Gaussian character of their marginal distribu- 
tions and negative values of the Wigner function. Other 
non-classical states have been recently produced starting 
from a squeezed vacuum and with the controlled subtrac- 
tion of a single photon: in this case the marginal distribu- 
tions clearly showed a squeezed and non-Gaussian char- 
acter but the preparation and detection efficiency was 
not high enough to reconstruct a negative- valued Wigner 
function 4]. 

We have recently reported the experimental gen- 
eration of a new kind of non-classical field states and 
their tomographic analysis based on time-resolved homo- 



dyne detection, which has allowed us to observe both the 
squeezed character and the negativity of the associated 
Wigner function. These so-called single-photon-added 
coherent states are produced whenever a single photon is 
injected in the same spatio-temporal mode of a coherent 
state and are shown to exhibit a mix of the character- 
istics of both parents. In particular, by simply varying 
the contribution of the initial coherent state, the charac- 
ter of the final state can be continuously tuned between 
completely quantum and almost completely classical. 

Here we report an in-depth tomographic analysis of 
such states based on new and more accurate experimen- 
tal data and comprising a complete comparison of the re- 
constructed density matrix elements with those expected 
from a fully developed theory. This has allowed us to fol- 
low in a much more detailed way the evolution of the gen- 
erated state from the particle-like one, characterized by 
a circularly symmetric and negative- valued single-photon 
Wigner function, through a squeezed intermediate region 
characterized by the gradual appearance of a phase, to- 
wards the wave-like classical coherent state. 

In addition to the interesting physical properties of 
single-photon-added coherent states, the ability to gen- 
erate, manipulate and characterize such states can be 
useful for possible future applications in the engineering 
of quantum states 0,0 and in quantum information pro- 
tocols _4J. 



PROPERTIES OF THE SPACS 

In 1991, Agarwal and Tara introduced a new class 
of states, defined by the repeated (m times) apphcation 
of the photon creation operator to the coherent state. 



— h 



-Am 



(1) 



2 



with ka^m = [to! Lm(— IcKp)]""^/^ & normalization fac- 
tor where Lm{x) is the mth-order Laguerre polynomial 
and m is an integer. Such photon-added coherent states 
(PACSs) essentially represent the result of successive el- 
ementary one-photon excitations of a classical coherent 
field and occupy an intermediate position between the 
Fock and the coherent states, reducing to the two limit 
cases for a — > or m — *■ 0, respectively. From the expan- 
sion of PACSs in terms of Fock states, it can be easily 
seen that they essentially correspond to a shifted version 
of a coherent state where all the \n) terms with n < m 
are missing and that all the elements of the correspond- 
ing density matrix are re-scaled and displaced towards 
higher indices pij — > pij^m.j+rm leaving all the elements 
with i,j<m void. 

When just a single quantum of field excitation is added 
to a coherent field, the singlc-photon-added coherent 
state (SPACS) reads as: 



\a) 



a- 



(2) 



and can be also viewed as the superposition of a dis- 
placed single-photon Fock state and a coherent state Q . 
SPACSs can be expanded in terms of Fock states as: 



Ia,l) 



ry'" 

^ V^min+l) (3) 



where the lack of the vacuum term contribution is evi- 
dent. Accordingly, the density matrix elements for the 
SPACSs are: 



Pi,j 



1) 



(4) 



and the effect of single-photon excitation can be readily 
observed in the plots of Fig.nwhere the absolute value of 
the theoretical matrix elements and the photon number 
distributions (their diagonal elements) are reported for 
the single-photon Fock state, for a coherent state with 
a = 1 and for the corresponding SPACS. 

Unlike the operation of photon annihilation, which 
maps a coherent state into another coherent state 
(d |a) — a ja)), i.e. a classical field into another classi- 
cal field, the single-photon excitation of a coherent state 
changes it into something quite different, especially for 
low values of a, where the absence of the vacuum term 
has a stronger impact. In the extreme case of an ini- 
tial vacuum state |0), the addition of one photon indeed 
transforms it into the very non-classical single-photon 
Fock state which exhibits negative values of the 
Wigner function around the origin. More generally, the 
Wigner function for a single-photon-added coherent state 
of arbitrary amplitude a can be expressed as: 



(5) 



a) 



b) 
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FIG. 1; (color online) Theoretical density matrix elements 
and photon number distributions for: a) the single-photon 
Fock state b) the coherent state \a) (with \a\ = 1) and 
c) the corresponding SPACS |a, 1) obtained by emission of a 
single photon in the mode of the coherent state. 



(where z = x ~\- iy) and can clearly become negative, a 
proof of its non-classical character, whenever the condi- 
tion 



\2z-aV < 1 



(6) 



is satisfied. Thus, in general, the application of the cre- 
ation operator d^, changes a completely classical coher- 
ent state into a quantum state with a varying degree of 
non-classicality, which becomes more evident the smaller 
the initial amplitude of the \a) state. If the amplitude 
a is gradually increased from zero, the smooth transi- 
tion from an initial purely quantum state (the single- 
photon Fock state) towards a classical coherent one (with 
the birth and the gradual appearance of a well defined 
phase) can be achieved. In addition to the negativity 
of the Wigner function, SPACSs also exhibit a definite 
squeezing in their field quadratures that can be readily 
observed. Given a field quadrature i:^ = ^(de~'^-|-d^e*^), 
its mean value is: 



(a, 11 Xe la, 1) 



|a|(2 + |a|2)cos(0) 



and its fluctuations amount to: 



1 1 



[Axo]l^{xl)^-{xe)i^- + 



lap cos (20) 



2(H-|a|2)2 



(7) 



(8) 



Clearly, the quadrature obtained by choosing 6 — ex- 
hibits reduced fluctuations with respect to the coherent 
state for |a| > 1, and is thus squeezed. It is interest- 
ing to note that, differently from Fock and Gaussian 
squeezed states, SPACSs combine both the key features 
normally associated to quantum states: the negativity of 
the Wigner function and the reduced fluctuations along 
one quadrature. 
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FIG. 2: (color online) Theoretical Wigner function for: a) the 
single-photon Fock state |1), b) the coherent state |a), c) the 
SPACS |q, 1). A value of \a\^ = 1 is used. 



GENERATION OF THE SPACS 

SPACSs can be generated by injecting a coherent state 
I a) into the signal mode of an optical parametric ampli- 
fier and exploiting the stimulated emission of a single 
down-converted photon into the same mode. Differently 
from conventional optical amplification where a coher- 
ent state is converted to another coherent state, here 
a low-gain regime of the amplifier and a conditioning 
of the state based on measurements on the idler mode 
are required to exactly select the one-photon excitation 
term (i.e. to avoid higher-order excitations which can- 
not be discriminated by our single-photon detectors and 
to exclude the vacuum contribution). Hence, in order to 
make sure that single-photon emission has taken place in 
the signal channel, one can use a conditional preparation 
technique which guarantees the generation of the target 
state every time that a single photon is detected in the 



correlated idler mode. 

The Hamiltonian for the parametric amplifier reads as: 

H = ihx{a\a\ - cisCLi) (9) 

where di (d^) is the annihilation operator for the idler 
(signal) mode and x is proportional to the amplitude of 
the (classical) pump and to the second-order susceptibil- 
ity of the medium. The time evolution of an initial state 
1^(0)) is thus described by: 

\m) = e-'"^ \m) - 

^gxt(ata+-a,a.) 1^(0)) . (10) 

If the parametric gain is kept sufficiently low (g = xt <C 
1), which is always the case in our experimental situation, 
the final output state can be approximated as 

m)) = [l + 9{alal-asa,)]\m)- (H) 

By letting a seed coherent field \a)^ enter the parametric 
crystal in the signal mode, while vacuum (|0)j) enters in 
the idler channel, the final state becomes 

\^) = [l+g{alal-aA)]\a)jO)^ = 

= |a)jO),+gdt|«)jl)^ (12) 

and the output signal mode will mostly contain the origi- 
nal coherent state, except for the few cases when the state 
|1) ■ is detected in the idler output mode. These relatively 
rare detection events, which take place with a probability 
proportional to |gp(l-|-|ap), project the signal state onto 
the desired SPACS corresponding to the stimu- 

lated emission of one photon in the same mode of \a). 
Note that when the input state is of the form \0)^ |0) -, i.e. 
no seed coherent field is injected into the crystal, spon- 
taneous parametric down-conversion takes place starting 
from the input vacuum fields, and pairs of entangled sig- 
nal and idler photons with random (but mutually cor- 
related) phases are produced in the crystal in the state 
|l)s ^'^^^ ^ low probability proportional to Igp. In 
this case, the detection of a single photon in the idler 
mode projects the signal state onto a single-photon Fock 
state and, by following the evolution of the final quantum 
state while the amplitude a increases from zero, one can 
witness the gradual transition from the spontaneous to 
the stimulated regimes of light emission with the smooth 
transformation of a single photon (particle-like) state to- 
wards a coherent (wave-like) one. 

Quite interestingly, one can obtain an absolute cali- 
bration of the amplitude of the seed coherent field Iq;)^ 
injected in the SPDC signal mode by measuring the rate 
of counts in the idler channel and comparing them to the 
un-seeded case. As stated above, the ratio of such rates 
equals (1 -I- lap) and this is clearly due to the enhance- 
ment of emission probability characteristic of stimulated 
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emission in bosonic fields. The same scheme was origi- 
nally proposed by Klyshko |^ as a metrological tool for 
absolute radiance measurements 0, 0| . 

For low a values, one can truncate the above expres- 
sions to the first two terms of the coherent state expan- 
sion in the number state basis. In this case the final state 
becomes: 

IV-) «[l+.9(alal-a.a,)](|0), + «|l)J |0), = 
= (|0), + «|1)J |0),+5(|l), + ^/2a|2)J (13) 

and the conditioning thus reduces the signal state to the 
coherent superposition + V2a 12)^,. Note that such a 
coherent superposition of two number states possesses a 
well defined phase and is highly non-classical, completely 
missing the contribution of the vacuum. 

The same state as the one described by (fT^ and 
has been recently generated and used by Resch et al. 
to generate an arbitrary superposition of zero- and one- 
photon states. In that case, however, the conditioning 
was performed upon the detection of a single photon in 
the same mode \\) ^ of the input coherent state, hence the 
final state was completely different from the ones investi- 
gated here, and of the form (a |0) -+3 The injection 
of a single photon instead of a coherent state as a seed 
for conditional parametric amplification has also been in- 
vestigated in 01 and experimentally demonstrated [l3| . 
with the amplification to a |2)^ Fock state in the context 
of quantum cloning. 

EXPERIMENTAL SETUP 

The experimental apparatus used to generate and ana- 
lyze the SPACS is schematically drawn in Fig|3| A mode- 
locked Tirsapphire laser, emitting 1-2 ps long pulses at 
786 nm and at a repetition rate of 82 MHz is used as the 
primary source. The laser pulses are frequency doubled 
to 393 nm in a 13-mm long LBO crystal which thus pro- 
duces the pump pulses for parametric down-conversion in 
a 3-mm thick, type-I BBO crystal. The crystal is slightly 
tilted from the collinear configuration in order to obtain 
an exit cone beam with an angle of ~ 3° from which 
symmetric signal and idler modes are roughly selected 
by means of irises placed at about 70 cm from the crys- 
tal. 

In order to non-locally select a pure state on the signal 
channel, idler photons undergo narrow spatial and fre- 
quency filtering before detection; indeed the nonlocally- 
prepared signal state will only approach a pure state if 
the filter transmission function is much narrower than 
the momentum and spectral widths of the pump beam 
generating the SPDC pair HI E [13, El- The idler 
beam is thus passed through a pair of etalon interference 
filters which perform a narrow (50 GHz) spectral selec- 
tion and is then coupled into a single-spatial-mode fiber 




FIG. 3; (color online) Experimental apparatus: HT-BS high 
transmission beam-splitter, LBO lithium triborate crystal, BS 
and BS-H 50% beam-splitters, VF variable attenuation filter, 
BBO-I type-I /3-barium borate down-converter crystal, PZT 
piezoelectric transducer, B.H.D. balanced homodyne detec- 
tor, F spectral and spatial filters, SPCM single photon count- 
ing module, LO local oscillator. 



before impinging onto a single photon counting module 
(Perkin-Elmer SPCM AQR-14). 

The weak coherent state a is obtained by controlled at- 
tenuation (VF in the figure, composed of a polarizer and 
a half-wave plate) of a small portion of the laser emission 
which is fed into the signal mode of the parametric crys- 
tal and is then directed to a 50% beam-splitter (BS-H 
in figure). Here it is overlapped with a second (intense) 
coherent state (again obtained from a portion of the orig- 
inal laser pulses) which is spatially mode-matched to the 
conditionally-prepared SPACS by the insertion of appro- 
priate lens combinations (not shown in the figure) along 
its path and serves as the local oscillator (LO) for the 
homodyne measurements [l9j . In order to finely adjust 
the alignment and the synchronization between the sig- 
nal and LO pulses, we use the stimulated beam produced 
by injecting a different seed pulse into the idler chan- 
nel of the parametric crystal. Under appropriate con- 
ditions [1^ . the beam generated by stimulated emission 
on the signal channel is emitted in a spatial mode which 
closely matches that of the target signal beam and can 
thus be used for alignment purposes. Measurements are 
performed at different values of the coherent seed ampli- 
tude |a| by rotating the half- wave plate; as seen above, a 
calibration of such an amplitude is simply obtained from 
a measurement of the increase in the idler count rate. 
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TIME-DOMAIN HOMODYNE MEASUREMENTS 

The pulsed homodyne detection scheme used to ana- 
lyze the quantum states has been recently developed by 
our group and is currently the only system capable of op- 
erating at the full repetition rate (80 MHz) of common 
mode- locked lasers in the time domain [3, [13 ■ The fields 
at the two output ports of the beam-splitter are detected 
by two photodiodes (Hamamatsu S3883, with active area 
1.7 mm^) whose difference signal is amplified and sent to 
a fast digital oscilloscope whose acquisition is triggered 
by the detection events in the idler channel. Each ac- 
quisition frame spans two consecutive LO pulses where 
only the first one is synchronized with the detection of 
an idler photon and contains the "information" about the 
SPACS \a, 1)^, while the second one can be used for the 
measurement of the reference un-excited coherent state 
|a)g. By blocking the seed coherent pulse, the single- 
photon Fock state [1)^, and the LO shot-noise distribu- 
tions corresponding to the vacuum state \0)^ are simul- 
taneously measured. About 5000 acquisition frames can 
be stored sequentially in the scope at a maximum rate of 
160,000 frames per second. Each sequence of frames is 
then transferred to a personal computer where the areas 
of the pulses are measured and their statistic distribu- 
tions are analyzed in real time. 

If a narrow temporal gate with the laser pulses is used, 
the typical rate of state preparation for vacuum input 
is about 300 s^^, with less than 1% contribution from 
accidental counts. A typical sequence of about 5000 ac- 
quisition frames can thus be captured and analyzed in 
about 20-30 s when no coherent seed is fed in the signal 
mode. It is interesting to remind that the probability 
of detecting an idler photon is proportional to |a^^ ja) p, 
hence, as soon as a is increased and stimulated emis- 
sion starts taking place, the rate of trigger events grows 
proportional to (1 + thus making the acquisition 

rate much higher. However, even at the maximum val- 
ues (I a I « 7) reached in the experiments, the trigger rate 
never exceeds 2 • 10'*s~^ (to be compared with the laser 
pulse repetition rate of about 8 • 10''s~^), so that the 
probability of conditioning the measurement upon more 
than a single idler photon always remains negligible. 

To explore the different quadratures of the generated 
field, its phase 9 relative to the LO field has been varied 
in controlled steps by applying a voltage to a piezoelec- 
tric transducer (PZT) which slightly translates one of 
the steering SPACS mirrors. A very good passive stabi- 
lization of the large (about 2 m long) Mach-Zehnder-like 
interferometer formed by the paths of the LO and of the 
seed coherent beam has been achieved in order to guar- 
antee a constant relative phase during the acquisition of 
a frame sequence. This is extremely tricky especially 
for low values of a, where the acquisition rate is lowest. 
About 50 acquisitions were performed at each value of 



the coherent seed amplitude a for 10-15 phase values in 
the [0,7r] interval. 

In order to reduce the contribution from low-frequency 
noise in the detection system, the amplified difference 
signal from the two photodiodes is AC-coupled (cut-off 
frequency of about 3 kHz at -3 dB) before subsequent 
amplification and acquisition. This high-pass filter, com- 
bined with the temporal sampling operated by the soft- 
ware for the measurement of the pulse areas, results in a 
reduced contribution from the DC terms in the acquired 
homodyne data which has to be taken into account in 
the analysis. While the fast pulse-to-pulse fluctuations 
which contribute to the marginal distributions are not 
affected by such filtering, the mean value of such distri- 
butions has to be scaled by a factor that depends on the 
sampling window used to measure the pulse area. This 
is clearly not an issue when measuring the vacuum field 
or Fock states, since the mean value of their marginal 
distributions is constantly zero, but it has to be carefully 
considered when dealing with states having a non-null 
mean field value. A calibration of such a factor can be 
simply carried out by comparing the sequence of marginal 
distributions as measured for the coherent state and the 
value of a as obtained from the idler counts, and can then 
be used to re-scale all the mean marginal values before 
subsequent analysis. This procedure was double-checked 
also by fitting the experimental SPACS marginals to the 
expected theoretical shapes and obtaining an indepen- 
dent measure of |a| to compare with the one deduced 
from the idler counts. 

DATA ANALYSIS AND DISCUSSION 

Balanced homodyne detection allows the measurement 
of the signal electric field quadratures xg = xcosO + 
y sin as a function of the relative phase 9 imposed be- 
tween the LO and the signal, where the two orthogonal 
field quadratures a; and y are defined as a; — ^{a+a^) and 
y = ^(a^ — d) and [x, y] — i/2. By performing a series of 
homodyne measurements on equally-prepared states it is 
possible to obtain the probability distributions p{x, 9) of 
the quadrature operator xg = i(ae~*^ -I- d^e^^) that are 
simply seen to correspond to the marginals of the Wigner 
quasi-probability distribution W{x,y) [2ll |: 

/ + 00 
W{x cos 9 — ysm9,xsin9 + y cos 9)dy. 
-oo 

(14) 

Given a sufficient number of quadrature distributions at 
different values of the phase 9 G [0, tt], one is able to re- 
construct the quantum state of the field under study . 
The elements of the density matrix p of the state in the 
number-state representation can be obtained by averag- 
ing the so called "pattern functions" fnm{x,9) over the 
outcomes of the quadrature operator and over the phase 
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8 as 



{n\p\m) = - 

IT 



+00 



dxp{x,0)fnm{x,9), (15) 



where the pattern functions can be implemented for unit 
quantum efhciency with stable numerical algorithms [23 . 
1231 . The Wigner function can then be obtained by means 
of the following transformation: 



W{x,y) 



M 

E 



Pn.rn Wn,m{x, y) 



(16) 



where Wn,m{x,y) is the Wigner function of the opera- 
tor \n) (m|. Note that, using this procedure, the Wigner 
function of the state is reconstructed from a truncated 
density matrix of dimension M x M . This implies a fi- 
nite resolution in the reconstructed function which, how- 
ever, can be adapted to the particular physical situation 
of interest in order to avoid loss of information on the 
state. 

In Fig. 0] a sequence of SPACS reconstructed density 
matrices is shown for increasing values of the seed co- 
herent field amplitudes \a\. In Fig. Othe corresponding 




FIG. 4: (color online) Density matrices of the SPACSs as re- 
constructed from the experimental data for increasing values 
of the seed amplitude, a) \a\ — i.e. single photon Fock 
state, b) |q| = 0.387, c) \a\ = 0.955, d) \a\ = 2.61. 

Wigner functions, obtained from the truncated density 
matrices of Fig 0] are shown: the first one (a), obtained 
with a blocked input, corresponds to the single-photon 
Fock state obtained by conditional preparation from the 
two-photon wavefunction of SPDC [ij, |2| and clearly ex- 
hibits classically impossible negative values around the 
center of the circularly symmetric (due to the undefined 
value of the phase) distribution. When the coherent seed 
is initially switched on at very low intensity (|a| « 0.4, 
i.e. an average of one photon every 7 pulses), the Wigner 
function starts to loose its circular symmetry while mov- 
ing away from the origin due to the gradual appearance of 
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FIG. 5: (color online) Wigner functions of the SPACSs as re- 
constructed from the density matrix elements shown in Fig0] 
Also shown are sections of the reconstructed Wigner functions 
in the a; = plane (data points) , together with the ones (solid 
lines) calculated as explained below by taking the limited ef- 
ficiency of the system into account. 



a defined phase, but it still exhibits a clear non-classical 
nature as indicated by its partial negativity (b). For 
increasing seed amplitudes, the negativity gradually gets 
less evident (c) and the ring-like wings in the distribution 
start to disappear making it more and more similar to the 
Gaussian typical of a classical coherent field (d). Inter- 
estingly, even at relatively high input amplitude a, the 
Wigner distribution for the SPACS \a, 1) keeps showing 
the effect of the one-photon excitation when compared 
to the corresponding, slightly displaced, un-excited la) 
state 0. 

When comparing the reconstructed Wigner functions 
and density matrix elements to the theoretical ones for 
the corresponding quantum states, one has to take into 
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account the limited efficiency of the homodyne detection 
apparatus which does not aUow one to generate and an- 
alyze pure states but always involves some mixing with 
the vacuum. The limited efficiency enters both in the 
preparation of the quantum state, where the dark counts 
and the non-ideal conditioning in the idler channel do 
not allow one to generate a completely pure state in the 
signal channel, both in the homodyne detection process 
itself, due to the limited efficiency of the photodiodes and 
to the imperfect mode-matching of the signal field with 
the LO 0. 

Here the limited efficiency can be measured directly by 
studying the single-photon Fock state obtained by block- 
ing the seed coherent field. The non-unit efficiency of 
the apparatus prevents the observation of the real single- 
photon Wigner function, and what one gets instead is its 
convolution with the vacuum one. The convolution re- 
sult is the well known s-parametrized quasi-probability 
distribution with the s parameter scaled by the detec- 
tion efficiency 77 |3| • From a fit of the experimental 
quadrature distributions to the corresponding theoretical 
phase-independent, marginal curves, we obtain an overall 
efficiency of 77 = 0.602 ± 0.002. 

The Wigner function of SPACSs in presence of limited 
efficiency is: 

121 



W{z) 



(17) 

and one can easily see that, as for single-photon Fock 
states, negative values can only be achieved with 77 > 0.5. 
The non-unit detection efficiency thus reduces the non- 
classical character of experimentally observed SPACSs 
and, especially for higher values of the seed amplitude \a\, 
may completely mask it in the presence of reconstruction 
noise (see Fig|SJi)). In our case, a good efficiency com- 
bined with relatively low reconstruction errors allow us 
to clearly observe the non-classical character of SPACSs 
up to \a\ « 2. 

The marginal distributions of the SPACSs in the case 
of limited efficiency have the form: 



1 



1 — ?7 -|- Arjx 



+ \a\^{l + 2t]{t] - 1)) - 4\a\xy/lj{2rj - 1) cos(6') 



+ 2|ap77(77 - 1) cos(26i)l e-2(:^-|"lv^cose) 



(18) 



and the corresponding mean values and variances have 
to be modified with respect to the ideal cases of Eqs.Q 
and ^ as: 



and 



|a|(2+ |Qp)^cos(6>) 

TTRp 

1 77(1 - |a|2cos(26')] 
4 + 2(l-f |a|2)2 



(19) 



(20) 



Figure|Blpresents the experimental marginal distributions 
of the SPACS at a fixed value of the coherent seed am- 
plitude \a\ = 0.387 and for three different values of the 
phase 9. Superposed to the data points are the curves 
obtained from a fit of the distributions to the expected 
shapes as given by ea. (|18|l . 
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FIG. 6: (color online) Normalized histograms of the pulse- 
integrated homodyne signal for the SPACS (filled circles) and 
the coherent seed field (empty circles) at different phases rel- 
ative to the LO. The quadrature x axis is normalized to the 
vacuum/coherent state distribution width. Also plotted (solid 
curves) are the fits to the theoretical curves for the Wigner 
marginals including the effects of the limited efficiency. 

The non-classical character of SPACSs is also evident if 
the quadrature variances are measured for different am- 
plitudes of the coherent seed field. Indeed, while the 
original coherent state has equal ffuctuations in the dif- 
ferent quadratures independently from its amplitude, the 
one-photon-excited state exhibits a squeezing in one of 
the quadratures and larger fluctuations in the orthogo- 
nal one as soon as \a\ > 1, as indicated in eqs.® and 
H20|l . An intuitive interpretation of this behavior can be 
connected with the reduction in the intensity noise of 
the coherent state when excited by a perfectly defined 
number of quanta with the corresponding increase in the 
phase noise due to the intrinsic lack of phase informa- 
tion of the Fock state. This effect starts to become ev- 
ident in the reconstructed Wigner function of Fig. El^c) 
(which is however still at the border of the un-squeezed 
region), where a somewhat reduced width appears along 
the radial direction, while the increase in the phase noise 
is indicated by the appearance of the ring-like wings in 
the tangential direction of the Wigner distribution. A 
more quantitative measurement of the variance in the 
squeezed and anti-squeezed quadratures is presented in 
Fig. [7| where the expected curves for = and 9 = Tr/2 
are also drawn according to eg. 1)20(1 with a global effi- 
ciency of 7; = 0.6. The experimental variances for the 
X(g=o) quadrature clearly get smaller than those of the 
corresponding coherent state (also shown in the graph 
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FIG. 7: (color online) Variances of the squeezed (filled 
squares) and anti-squeezed (filled circles) quadratures of the 
SPACS for different coherent state amplitudes. Solid and 
dashed lines are obtained from Eq. 1201 with 6 — and 
9 — Tr/2, respectively, and with a global efficiency set to 
77 — 0.6. Also shown are the experimental data (empty cir- 
cles) and the theoretical curve (horizontal line at 1/4) for the 
variance of the coherent state. 



and independent of the seed intensity) as soon as the 
amplitude exceeds unity, and a maximum squeezing of 
about 15% is obtained for \a\ — 1.85. 

The density matrix elements reconstructed from the 
data can also be compared with the theoretical ones pro- 
vided that a Bernoulli transformation 



fc=0 



J + k 



-.1/2 



(l-ry)V.,, (21) 



is performed in order to include the effects of non-unit ef- 
ficiency. The expected density matrix pc is thus obtained 
from Eqs.lQ and with rj = 0.6. Figure |H1 shows the 
calculated density matrix elements and the correspond- 
ing Wigner function for the SPACS with \a\ — 0.955. 
Such data should be compared to the experimental plots 




FIG. 8: (color online) a) Calculated density matrix elements 
and b) Wigner function of the SPACS in the case of limited 
efficiency (77 = 0.6) and for \a\ — 0.955. 

of Figs. Et) and Et). 



The theoretical and experimental photon number dis- 
tributions for the SPACS with increasing seed amplitudes 
\a\ are also plotted in Fig. El where the errors have been 
calculated as in Ref.'23l. 
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FIG. 9: (color online) Experimental (data points with error 
bars) and theoretical (solid bars) photon number distributions 
for the SPACS in the case of limited efficiency (77 = 0.6). a) 
|q| = 0, b) \a\ = 0.387, c) \a\ = 0.955, d) \a\ = 2.61. 

In order to better quantify the agreement of the exper- 
imental results with the theoretical ones, we calculated 
the purity P of the reconstructed states and compared it 
with the expected one Pc- The purity tends to increase 
with the amplitude of the seed coherent field as the fi- 
nal state evolves from a mixture of vacuum and single 
photon towards a pure coherent one. Results are pre- 
sented in Table HI where we also show the fidelity of the 
reconstructed states to the theoretical ones, calculated as 
proposed by Jozsa |25j for the comparison between mixed 
quantum states 



F = 



Tr 



IPcPe\JPc 



(22) 



where pe is the experimentally derived density matrix. 



TABLE I: Quantum state indicators for SPACSs. M dimen- 
sion of the reconstructed density matrix; Pc expected purity 
for the state; P purity of the reconstructed state; F fidelity 
between the experimentally reconstructed and the expected 
state. 



a 


M 


Pc = Tr(p2) 


P = Tr{pl) 


F 





6 


0.52 


0.53±0.01 


1.009±0.004 


0.387 


7 


0.64 


0.63±0.01 


0.990±0.004 


0.955 


8 


0.87 


0.83±0.02 


0.976±0.004 


2.61 


14 


0.99 


1.02±0.05 


0.995±0.005 



A comparison of the fidelity values with the photon 
number distributions of Fig. El seems to indicate that fi- 
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delity, quickly saturating to unity even for somewhat dif- 
ferent distributions, is not a very sensitive benchmark for 
the closeness of the experimental to the expected states 
in this case. A similar conclusion has been recently ob- 
tained by the group of Kwiat |26j in the case of depo- 
larized entangled mixed states in the discrete variable 
domain. 

CONCLUSIONS 

We have successfully used a conditional preparation 
technique to generate a new class of light states whose 
degree of non-classicality can be continuously tuned be- 
tween the extreme situations of pure quantum states 
and classical ones. Such singlc-photon-added coherent 
states are particularly interesting from a fundamental 
point of view as they represent the result of the most 
elementary excitation of a classical light field. In this 
regard, the demonstrated possibility to follow their evo- 
lution so closely will certainly push the experimental re- 
search towards the investigation of other interesting and 
equally fundamental quantum processes. Single-photon- 
added coherent states are also noteworthy because, for 
the first time to our knowledge, both the typical prop- 
erties of a quantum character, i.e. the negativity of the 
Wigner function combined with a quadrature squeezing, 
are simultaneously and clearly observed in a light state. 
While the non-classicality criterion based on the squeez- 
ing can be in principle fulfilled without bounds on the 
detection efficiency, as shown in ea.()20|l. the negativity 
of the Wigner function is a much stricter criterion since 
it requires efficiencies higher than 50%. We have demon- 
strated the possibility of using our recently developed 
technique for high-frequency time-resolved balanced ho- 
modyne detection to reconstruct the density matrix ele- 
ments and the Wigner functions of the generated states 
with an overall efhciency of 60%. Thanks to the very 
high acquisition rates achievable by our system, other 
non-classical states, particularly those involving higher 
number of photons and normally characterized by lower 
generation efficiencies, are within reach for a complete 
tomographic analysis. 
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